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ABSTRACT 

Type Il-plateau supernovae (SNe IIP) are the results of the explosions of red supergiants and are 
the most common subclass of core-collapse supernovae. Past observations have shown that the outer 
layers of the ejecta of SNe IIP are largely spherical, but the degree of asphericity increases toward 
the core. We present evidence for high degrees of asphericity in the inner cores of three recent 
SNe IIP (SNe 2006my, 2006ov, and 2007aa), as revealed by late-time optical spectropolarimetry. The 
three objects were all selected to have very low interstellar polarization (ISP), which minimizes the 
uncertainties in ISP removal and allows us to use the continuum polarization as a tracer of asphericity. 
The three objects have intrinsic continuum polarizations in the range of 0.83 — 1.56% in observations 
taken after the end of the photometric plateau, with the polarization dropping to almost zero at the 
wavelengths of strong emission lines. Our observations of SN 2007aa at earlier times, taken on the 
photometric plateau, show contrastingly smaller continuum polarizations (~ 0.1%). The late-time 
Ha and [O I] line profiles of SN 2006ov provide further evidence for asphericities in the inner ejecta. 
Such high core polarizations in very ordinary core-collapse supernovae provide further evidence that 
essentially all core-collapse supernova explosions are highly aspherical, even if the outer parts of the 
ejecta show only small deviations from spherical symmetry. 

Subject headings: polarization — supernovae: individual (SN 2006my, SN 2006ov, SN 2007aa) 



1. INTRODUCTION 

Core-collapse supernovae (SNe) are the explosive 
deaths of massive stars. Indirect inferences that these 
energetic explosions are aspherical have come from such 
evidence as the hi gh space velocities of t he pulsars pro- 
duced by the SNe ([Lvne fe Lorimerlll994D and the struc- 
ture of Galactic supernova remnants (e.g., Fesen 2001). 
A direct observational tes t for the presence of asphericit y 
in SNe was proposed by IShapiro fc Sutherland! (|1982f ). 
who pointed out that the polarizing effects of electron 
scattering in a supernova atmosphere would not cancel 
if the supernova were aspherical and spatially unresolved 
by a distant observer. 

The first published dete ctions of supernova polariza- 
tion were for SN 1987A (jBarrettl 119881 : ICropper et all 
119881: iMendez et aTlll988[ ). Subsequent observations of 
additional SNe led to the discovery that all core-collapse 
supernovae with data of suffici ent quality exhibit polar - 
ization, as recently reviewed by I Wang fc Wheeler! (2008). 
An important test of our understanding was provided 
by the late-time Hubble Space Telescope (HST) images 
of SN 1987A which spatially resolve the ejecta and 
show that they are elongated along a position angle 
(P.A.) consistent w ith that inferred from the polarimetry 
([Wang et alJ(2002]) . 

Our study here focuses on SNe IIP, which represent 
the single most common outcome for the ends of the 
lives of massive stars. Approximately 50% of all core- 
collapse SNe are of this subtype (Li et al., in prep.; see 
also Smartt et al. 2009). Their progenitors are rel- 
atively low-mass red supergiants in the mass range of 
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8.5-16.5 M Q ([Smartt et alJ"2 009) that are either isolated 
or sufficiently separated from their binary companions so 
as to leave their evolution and outer envelopes unaffe cted 
([Heger et alJl2003tlLi et"aT]|2006l:[slnartt et al.ll2009h . In 
many ways, their progenitors and explosions should be 
the simplest ones to understand and model. 

Despite this, published spectropolarimetry of SNe IIP 
is surprisingly limited. Mos t of the objects of this t ype 
listed in the compilation of Wa ng fc Whe eler (200||) as 
having available polarimetry are either unpublished or 
have only photometric polarimetry. Our focus here is on 
measuring continuum polarization and separating intrin- 
sic continuum polarization from the effects of interstellar 
polarization (ISP); with only photometric data, this is 
extremely difficult, if not impossible. 

The two SNe IIP that have been the most well-studied 
spectropolarimetrically are SNe 1999em and 2004dj. 
SN 1999em showed low polarization shortly after explo- 
sion, indicating that the outer ejecta were nearly spher- 
ical. The level of continuum polarization slowly in- 
creased with time in a manner consistent with a sin- 
gle axis of symmetry while simulta neously the po l ariza- 
tion signatures in the lines grew ([Leonard et al.l 120011 : 
IWang fc Wheelerll2008t) . SN 2004dj showed no detectable 
polarization (after correction for ISP) in the continuum 
or lines while on the plateau. At the end of the plateau 
stage, the continuum polariz ation jumped up to 0.56%, 
revealing an aspherical core (|Leonard et~aI1l2006h . The 
polarization percentage slowly decreased thereafter in a 
manner attributed to geometric dilution of the scatter- 
ing electrons as the ejecta expanded. The observations of 
SN 2004dj prompted the present study in an attempt to 
determine whether the results for SN 2004dj were com- 
mon to SNe IIP. 

Three recent SNe IIP in the Virgo Cluster have pro- 
vided an opportunity to study this class of objects in 
more detail. SN 2006my was discovered by K. Itagaki 
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in NGC 4651 on 2006 November 8.82 (UT dates are 
used throughout this paper; Nakano & Itagaki 2006a). 
The same observer also discovered SN 2006ov in M61 on 
2006 November 24.86 (jNakano fc Itagakil |2006bl) . Both 
objects were soon spectroscopically confirmed as SNe IIP 
a few months past explosion (Stanishev & Nielsen 2006; 
Blondin et al. 2006). SN 2007aa in NGC 4030 was dis- 
covered by T. Doi on 2007 February 18.308 and spec- 
tra obtained by N. Morrell on Feb. 19.24 showed that 
it was a SN IIP simil ar to SN 1999em at about 2 
days after expl osion (|Doll2007t iFolatelli fe Morrellll2007D . 
iLi et al.l (|2007ft presented photometric and spectroscopic 
evidence that both SN 2006my and SN 2006ov were nor- 
mal Type IIP SNe discovered near the end of the photo- 
metric p lateau st a ge. 

While ILi et al.l (|2007|) identified possible progenitors 
for SNe 2006my and 2006ov in pre-cxplosion HST im- 
ages, subsequent work has rejected thos e identifications 
dLeonard et al.ll2009t ISmartt et all l2009t [Crockett et al 
12009ft Instead. ISmartt et all (|2009h and ICrockett etTal 
(2009) give upper limits for the zero-age main-sequence 
masses of any progenitor in the HST images of 13, 10, 
and 12 M q for SNe 20 6my 2006ov, and 2007aa, re- 
spectively. ISmartt et alJ (|2009| ) also quote values for the 
56 Ni masses produced by SN 2006ov and SN 2006my 
of 0.003 ± 0.002 Mq and 0.03 ± 0.015 M Q , respectively, 
based on analyses of the late-time light curves. The 56 Ni 
mass inferred for SN 2006ov is extremely low, compa- 
rable to that of the underlu minous SN 1999br (|Hamuvl 
120031: iPastorello et~all l2004l ) . despite SN 2006ov having 
a normal luminosity on the plateau. 

In this paper, we present spectropolarimetry of the 
above three SNe IIP obtained after the end of the photo- 
metric plateau. In §2, we describe the observations and 
corrections for small amounts of ISP from the Galaxy. 
We present the evidence for very high intrinsic contin- 
uum polarizations in these objects in §3 and discuss the 
interpretation and implications of our results in §4. 

2. OBSERVATIONS 

Unfiltcrcd images of all three SNe in our sample were 
obtained as part of our regular monitoring of their host 
galaxies for the Lick Observatory Supernova Search us- 
ing the 0.76-m Katzman Automatic Imaging Telescope 
(KAIT; Filippenko et al. 2001), although we artificially 
increased the frequency of observations from the default 
search strategy to better sample the light curves of these 
objects. Portions of the SN 2006ov and SN 2006my light 
curves have a lready b e en pre sented and the reductions 
described by ILi et alJ (|2007l ). Briefly, a point-spread 
function fitting technique was used to subtract pre-SN 
template images from the SN observations and aper- 
ture photometry was p erformed on t he subtracted im- 
ages using DAOPHOT (jStetsonl Il987t ) in IRAF 3 . KAIT 
unfiltered obse rvations approximate the standard R band 
(|Li et al.ll2003l ). so the SN observations were then cali- 
brated to an average of USNO-B1 stars in the field (up to 
20 per field), resulting in systematic uncertainties in the 
zero points for each field of 0.05 — 0.18 mag. However, 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories (NOAO), which are operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



only the shapes of the light curves are relevant for this 
study. 

We obtained spectopolarimctry of the three SNe us- 
ing both t he Kast spectrograp h on the Lick 3-m Shane 
telescope ( Miller fe Stond[l993| ) and the Low Resolution 
Imaging Spectr ometer (LRIS) m ounted on the Keck I 
10-m telescope (|Oke et al.lll995f) . The Kast data were 
obtained using the red camera only, with the 300/7500 
grating tilted to cover a wavelength range of approxi- 
mately 4400—9900 A in conjunction with a GG455 order- 
blocking filter. A small amount of second-order light is 
likely present at wavelengths greater than ~ 9000 A, 
but none of our conclusions depend on data at those 
wavelengths. A 3"-wide slit gave a spectral resolution 
of ~ 16 A. 

The LRIS data were taken with both halves of the 
spectrograph. All observations used an identical setup 
with the D560 dichroic beam-splitter to split the spec- 
trum near 5500 - 5600 A, along with the 600/5000 grism 
on the blue side and the 400/7500 grating on the red side. 
We used a l'.'5-wide slit to achieve spectral resolutions 
of ~ 6 A and 9 A on the blue and red sides, respec- 
tively. The total spectral coverage was 3170 — 9240 A, 
but the polarization data plotted in this paper are typi- 
cally truncated below ~ 4000 A because the low blue flux 
in SNe IIP due to line blanketing, and decreasing detec- 
tor sensitivity at short wavelengths, combine to make the 
polarimetry extremely noisy at the blue end. In addition, 
this spectral region suffers from background-subtraction 
uncertainties due to the underlying blue starlight from 
the host galaxy. Further details of the individual obser- 
vations are given in Tabic [T] 

Basic two-dimensional image processing and spectral 
extraction were accomplished using standard tasks in 
IRAF. Flux cal ibration and removal of atmospheric ab- 
sorption bands (|Wade fe Ho rnc 1988 ) was perf ormed us- 
ing our own IDL tasks (jMatheson et al.ll2000l ). The in- 
dividual one-dimensional SN spectra were combined on 
a pixel-by-pixel basis, weighted using the error est imates 
from the optimal spectral extraction technique (jHornd 
1986), and the final blue-side and red-side spectra (for 
LRIS) were scaled and combined across the overlap re- 
gion. In those cases where differential loss of blue light 
was apparent because the observed P. A. of the obser- 
vatio ns deviated from the parallactic angle (Filippenko 
1982), the combined observation was warped to the shape 
of the observation obtained closest to the parallactic an- 
gle using a low-order polynomial fit to the ratio of the 
spectra. 

We also extracted some of the emission from nearby 
H II regions along the slits of selected observations. We 
fit Gaussians to the nebular emission lines and obtained 
recession velocities of 920, 1660, and 1650 km s _1 for the 
H II regions near SN 2006my, SN 2006ov, and SN 2007aa, 
respectively. These values are within a typical value for 
a galactic rotation velocity (100 — 200 km s _1 ) of those 
listed for the host-galaxy nuclei in the NASA Extragalac- 
tic Database (NED), so we adopt these redshifts in this 
paper and have removed them from all plots. All refer- 
ences to wavelengths are in the rest frame of each object. 

2.1. Spectropolarimetry 

The spectropolarimetric obser vations were r e duced fol- 
lowing the standard method of iMiller et alJ ()1988l ). as 
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TABLE 1 
Details of Observations 



Object 


UT Date 


Age a 


Tclcscope b 


Total Exposure Time 


Seeing 


P.A. 


Notes 




(YYYY MM DD.DD) 


(d) 




to 


(") 


(°) 




SN 2006ov 


2006 12 20.51 


-5 


L3 


4800 


2.2 


144.5 


clouds 


SN 2006ov 


2006 12 25.63 





Kl 


4000 


1.4 


125 


clear 


SN 2006ov 


2007 01 21.56 


27 


Kl 


4400 


0.8 


125 


clear 


SN 2006my 


2007 01 21.50 


38 


Kl 


3600 


0.9 


282 


clear 


SN 2007aa 


2007 03 18.45 


-50 


Kl 


3600 


1.0 


359.7 


clear, poor guiding 


SN 2007aa 


2007 04 15.48 


-22 


Kl 


2800 


1.8 


57.5 


clear 


SN 2007aa 


2007 05 10.31 


3 


L3 


7200 


2.1 


203 


clear 



a Age relative to the end of the plateau, t p (see §3 for a definition). 
" L3 = Lick 3-m + Kast spectrograph; Kl = Keck I 10-m + LRISp. 



implemented by iLeonard et al.l (|2001l) . The angle zero 
points of the waveplat es were determined u sing polari- 
metric standards from ISchmidt et "all (|1992t ). Typically, 
two standards were observed per night to confirm the 
zero point. In addition, we noted that the zero points 
for each instrument were generally stable over time 
with small shifts. Low-polarization s t andards from the 
litera t ure (iMathewson fc Fordl Il970t lei cmens fc Tapial 
1990; Schmid t et al.lll992t iBerdvugin et al.l fl995) were 
observed on each night to have negligible polarizations, 
setting a limit on the instrumental polarization of P < 
0.1%. 

This study hinges on our ability to accurately measure 
the continuum polarization of SNe. ISP in both our own 
Galaxy and the host galaxies of these SNe can mimic 
the presence of continuum polarization and contaminate 
our results. All three SNe selected for this study are at 
high Galactic lat itudes (b > 59°) and h ence low Galactic 
ISP is expected (jSerkowski et al.lll975| ) based on the low 
columns (E( B — V) < 0.03 mag) measured along these 
lines of sight (|Schlegel et al.lll998l ). We can directly mea- 
sure the G alactic contrib ution to the ISP by following 
the lead of iTranl (|1995al ) and selecting stars of low in- 
trinsic polarizations and sufficiently large distances (es- 
timated via spectroscopic parallaxes) to place them out 
of the Galactic plane. Such probe stars should sample 
most of the Galactic contribution to the ISP along these 
lines of sight. The probe stars for SN 2006my were BD 
+ 17°2526 and BD +16°2409, those for SN 2006ov were 
BD +05°2618 and BD +04°2608, while BD +00°2876 
and BD -00°2514 were selected for SN 2007aa. 

The probe stars were observed immediately after the 
object on one observation date per SN with the slit kept 
at the same P.A. Each pair of stars gave consistent re- 
sults (within 0.1% in each Stokes parameter), so the two 
stars were averaged to form our estimate of the Galactic 
ISP for each SN. As expected, the Galactic ISPs were 
very small, with F-band polarizations of 0.07%, 0.06%, 
and 0.16% measured for the SN 2006my, SN 2006ov, and 
SN 2007aa probes, respectively. For the sake of com- 
pleteness, we subtracted low-order polynomial fits to the 
Stokes q and u of the combined probes from the Stokes 
parameters measured during each SN observation, al- 
though the effect is almost negligible for SN 2006my and 
SN 2006ov. 

As a final check on our instrumental stability, we 
observed the same probe star, BD +05°2618, on all 
three occasions that we obtained spectropolarimetry of 
SN 2006ov. The variation in each Stokes parameter 



was < 0.04%, consistent with similar experiments we 
have performed in the past with repeated observations 
of bright stars on different nights. This sets an upper 
limit to any systematic offset between the Lick and Keck 
polarimeters as well as to an y night-to-night instrumen- 
tal instability. We note that ILeonard et al.l (|2006l ) also 
found similar results using the Kast polarimeter. There- 
fore, we adopt a value of 0.04% for any systematic error 
contribution where appropriate in the analysis below. 

3. RESULTS 



As shown by iLi et all (|2007| ) for SN 2006my and 
SN 2006ov, and below for SN 2007aa, all three objects 
in our study exhibit standard light-curve shapes for SNe 
IIP, with the three phases of a relatively flat plateau 
followed by a rapid decline phase and then a steady de- 
cline at a rate consistent with the radioactive decay of 
56 Co. Thi s characteristic shape is well understood theo- 
retically (jEastman et al.l 11994 iKasen fc Wooslev l 12009 ) 
and defines this subcla ss of SNe (|Barbon et ah I197S ; 
iDoggett fc Branchlll985fh The supernova shock deposits 
thermal energy in the envelope which is released as the 
hydrogen in the envelope recombines during the plateau 
phase. When the photosphere recedes all the way back 
through the hydrogen layer, recombination ceases to be 
a source of energy and the SN light curve rapidly falls 
until the luminosity matches the instantaneous rate of 
energy input from the decay of 56 Co. 

Both SN 2006my and SN 2006ov were discovered 
several months after explosion near the ends of their 
plateaus. In order to make temporal comparisons 
amongst our observations of these SNe when the dates of 
explosion are so poorly known, we found it convenient to 
define a common time to represent the end of the plateau. 
As discussed below, the timing relative to the end of the 
plateau is probably the important physical parameter for 
t he polarization , not the time elapsed since explosion. 

Hamuy (2003) defined the end of the plateau phase, t p , 
as the time when the SN magnitude falls halfway between 
its value on the plateau and the value at the start of 
the radioactive decay tail. The plateau magnitude can 
be difficult to uniquely define in some cases due to the 
varied shapes of SN IIP light curves. Th e comp ilation 
of SN IIP photometry by iPoznanski et al.1 ([2009) shows 
that some objects have "plateau" phases that actually 
rise for many days before falling while others steadily, 
if slowly, decline from the beginning. The photometric 
band used to define the light curve can also have an effect 
as the shapes of SN IIP light curves are a function of 
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wavelength. iHamuvl (|2003l) used the V band, but we 
have only unfiltered photometry for our three objects. 

Despite these uncertainties in the definition of the 
plateau magnitude, t p occurs during the rapid decline 
phase and hence the uncertainty in t p is at worst a few 
days, which is negligible for our purposes. We measure 
dates for t p of 2006 Dec. 15.0, 2006 Dec. 25.5, and 2007 
May 7.3 for SNe 2006my, 2006ov, and 2007aa, respec- 
tively, based on the KAIT unfiltered photometry. Below, 
we measure t p for SN 1999em in the R band (which most 
closely approximates our KAIT unfiltered photometry), 
and our value f or t v differs b y only 0.2 days from the 
value derived by IHamuvl (|2003l ) in V, so we can conclude 
that bandpass effects have negligible effect on our mea- 
sured values for t p . Hereafter, dates listed as "day X" 
refer to epochs of observation X days after the end of the 
p lateau as define d by t p . 

iLi et al.1 (|200l found that SN 2006my and SN 2006ov 
were not highly reddened and had no detectable Na I D 
absorption lines from the interstellar media of their hosts. 
We confirm the lack of narrow Na I absorption in our 
spectra. Our highest signal-to-noise ratio (S/N) spectra 
of each object are plotted in Figure Q] SN 2007aa does 
have a weak dip near 5888 A that may be from absorption 
in NGC 4030. The strength of the feature in SN 2007aa 
is difficult to quantify at our low spectral resolution, es- 
pecially with the curvature of the spectrum making the 
continuum hard to define. We estimate the equivalent 
width of the dip to be ~ 0.08 A, with an uncertainty of 
about 50%. If two-thirds of the absorption is due to the 
Dl component of the doublet (reflecting the ratio of the 
oscillator strengths, as ex pected in the low optical depth 
limit), the relationship of iMunari fc Zwitterl lfl99l pre- 
dicts a host-galaxy reddening of E{B — V) « 0.014 mag, 
which corresponds to an ISP of P w 0.04% for typical 
lines of sight in the Galaxy, up to an e mpirical maximum 
of P « 0.13% (ISerkowski et al.l[l975l) . SNe 2006my and 
2006ov appear to have even less absorption from their 
host galaxies. Even if our estimates of the reddening are 
off by a factor of a few, the conclusion that we should 
expect our objects to have low host-galaxy ISP will be 
unchanged. 

3.1. SN 2006ov 

We obtained three epochs of spectropolarimetry on 
SN 2006ov, two during the steep fall off the photometric 
plateau and one a month later. Our highest S/N data 
were obtained from Keck on day and were fortuitously 
exactly synchronized with the end of the plateau. These 
data are plotted in the q ~ u plane in Figure [5J with 
each black circle representing a 100 A bin of data. As 
described above, the S/N degrades to the blue, so we 
only plot data at wavelengths greater than 4500 A. The 
cloud of points is clearly elongated in the positive q, pos- 
itive u direction by an amount in excess of the error bars, 
indicating substantial intrinsic supernova polarization. 

We selected two broad wavelength intervals, 6800-7200 
and 7820-8140 A, as representative of the line-free con- 
tinuum (the exact choice is justified below) and plotted 
the integrated polarizations 4 as the blue diamond and 

4 All integrated polarizations in this paper are calculated from 
the flux-weighted averages of q and u over the indicated wavelength 
intervals. 
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Fig. 1. — The spectral region near the Na I AA5890, 5896 dou- 
blet for each of the three SNe. The plotted spectra are from days 
38, 0, and -50 for SNe 2006my. 2006ov, and 2007aa, respectively. 
SN 2007aa shows a tiny dip near 5888 A that may be Na I absorp- 
tion in its host galaxy. The other two SNe show even less evidence 
for host-galaxy absorption, indicating that we should expect neg- 
ligible ISP from the host galaxies of these objects. 

solid red circle in Figure [5] The two continuum points 
are very consistent with each other and have an average 
polarization (weighted by the inverses of their variances) 
of P = 1.56 ± 0.03% at P.A.= 8?2 ± 0?6. Many of the 
100 A bins of data are clustered around these continuum 
points and along the line connecting them to the origin. 
The P. A. of 8? 2 therefore is significant in this object and 
probably represents the projection of the SN symmetry 
axis onto the plane of the sky. Since the Stokes param- 
eters (q,u) are referenced to the arbitrary direction of 
North on the sky, it is useful to rotate our coordinate 
system to align one Stokes parameter with the symme- 
try axis of the SN. We call these new coordinates rotated 
Stokes parameters (RSP; Trammell et al. 1993a; Tran 
1995b). If the SN were axisymmetric, ^rsp would be an 
estimator of the total polarization and ursp would show 
any deviations from axisymmetry. 

The polarization data after rotation are plotted in Fig- 
ure [3] in 20 A bins. The most striking aspect of Figure [3] 
is the high level of polarization (P rs 1.5%) in the con- 
tinuum, with many strong, sharp depolarizations down 
to almost zero in Qrsp present at the lines. We inte- 
grated the polarimetry over the wavelength range 6540- 
6600 A to determine the polarization at the peak of Ha 
and measured (grtsp, ursp) 01 (0.08%, 0.00%), with er- 
ror bars of 0.05% in each Stokes parameter. The similar 
depolarizations present for many other lines are convinc- 
ing evidence that the host-galaxy contribution to the ISP 
must be very small (< 0.1%), as might be expected from 
the low reddening to SN 2006ov. To estimate the ISP, 
many past workers have made the assumption that the 
intrinsic pola rization in supernovae at the peaks of strong 
lines is zero (iTrammell et al.lll993b1 : iHofhch et alJll996t 
iTran et aTlll997D. but that assumpt ion is clearly wrong 
in some cases (|Maund et aT1l2007bl ). especially at early 
times when He I A6678 is blended with Ha. Here, the 
continuum polarization we measure must be almost en- 
tirely intrinsic to SN 2006ov and the lines are in fact 
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Fig. 2. — SN 2006ov data from 2006 Dec. 25 (day 0) plotted in the q — u plane. The black circles represent 100 A bins of data at 
wavelengths greater than 4500 A. The radii of the circles are proportional to wavelength. The blue diamond and solid red circle represent 
the polarization in two continuum bins while the purple triangle marks the peak of the Ha emission line. The dashed lines represent the 
preferred, rotated set of Stokes parameters c/rsp and mrsp, chosen to align <jrsp with the two continuum bins. The cross at (—0.25,-1.3) 
shows the size of the median error bars on the data points. 



depolarized. 

With that in mind, we can now justify our choice of 
continuum intervals plotted in Figure [5] and marked on 
Figure [3l Past work has noted the utility of measur- 
ing the continuum polarization in the red portion of the 
optical w here few line features interfere with the mea- 
surement. iLeonard et al.l (|2006l ) used a wide continuum 
region of 6800-8200 A in their study of SN 2004dj. In- 
tegrating the observed polarization over a continuum re- 
gion that wide for SN 2006ov will clearly underestimate 
the true polarization as it will include a strong depo- 
larization present at the K I/O I blend as well as a 
weaker one due to [Ca II] emission. We formally get 
P = 1.26±0.02% if we use the wide definition for the con- 
tinuum, significantly less than the value of 1.56% quoted 
above. Therefore, we define two continuum regions to 
exclude the strong lines. The first one, 6800-7200 A, is 
designed to avoid the red wing of Ha at its blue end and 



also exclude the [Ca II] emission at the red end. The 
second region, 7820 - 8140 A, is bracketed by O I A7774 
and Na I AA8183, 8194. 

We also obtained SN 2006ov data on two additional 
epochs, days —5 and 27. We have plotted a comparison 
of <Zrsp on all three epochs in Figure |3J Due to the lower 
S /N of the additional data, we have rebinned to 50 A per 
pixel. Despite the significant noise, the overall pattern is 
the same: large continuum polarization with dips down 
to nearly zero present at the wavelengths of strong line 
features. The Lick data were taken only 5 days before 
the day Keck data and do show some significant dif- 
ferences, indicating rapid evolution of the spectrum as 
the SN made the transition off the photometric plateau. 
For example, the depolarizations in the wavelength range 
6000-6500 A are more blended together at day —5 than 
at day 0. 

We calculated the polarization in both of our contin- 
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Fig. 3. — Spectropolarimetry of SN 2006ov on day 0, presented as 20 A bins in the rotated (?r,sp~ "rsp coordinate system so that <zrsp is 
an estimator of the polarization and mrsp is near zero. The total-flux spectrum of SN 2006ov is overpolotted in the top panel to guide the 
eye, with several major lines labeled. The continuum polarization between strong line features is very high (~ 1.5%), but the lines result in 
depolarization. The two horizontal black bars mark the extent of our two relatively line-free regions (6800-7200 and 7820-8140 A), which 
we take to be representative of the continuum. 

uum regions for each of our three epochs. The results 
are plotted in Figure [5] and compared to the KAIT un- 
filtered light curve. The two continuum windows gave 
very consistent results, showing a significant increase in 
the weighted-average polarization in only five days, from 
1.25±0.06% on day -5 to 1.56±0.03% on day 0, falling 
to 1.15±0.08% on day 27. SN 2004dj showed large vari- 
ations in the polariza tion angle at the end of the plateau 
([Leonard et al.ll2006| ). but the SN 2006ov continuum ob- 
servations are consistent with a single P.A. to within 1°. 
Unfortunately, SN 2006ov exploded while in solar con- 
junction and we were unable to obtain observations in 
the month between discovery and our day —5 data, so 
it was impossible to establish whether the high polar- 
ization apparent at late times represented a significant 
change from early times. 

3.2. SN 2006my 

We obtained a single epoch of spectropolarimetry of 
SN 2006my on day 38, 2007 Jan. 21. We analyzed the 
data in a similar manner to SN 2006ov above. First, we 
measured the weighted average polarization in our two 
continuum windows to be P = 0.97 ± 0.04% at P.A. = 
174?1 ± 1?2. Then we rotated our Stokes parameters by 
174?1 to align <?rsp with this polarization angle. The 
results are plotted in Figure |6l 

The data are of lower quality than for SN 2006ov 
but show broadly similar characteristics. High contin- 
uum polarization is seen in (?rsp, with deep depolariza- 
tions present a t the wavelengths of strong line features. 
iLi et all ([20071 ) found that the extinction to SN 2006my 
was likely to be minimal, and hence low ISP should be 
expected. We determined the polarization at the peak 
of Ha using the same window as for SN 2006ov above 
and found (<7rsp,mrsp) of (0.03%, 0.03%) with an un- 
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Fig. 4. — Three epochs of polarization data for SN 2006ov. The 
data points with error bars represent the Stokes parameter QasPi 
which is an estimator of the total polarization. As a reference, 
the solid line in each figure is the total-flux spectrum from each 
epoch. Note the high continuum polarization, with depressions 
at the location of the strong lines. The high point near 8150 A 
in the day —5 polarization is spurious and due to relatively poor 
night-sky subtraction. The two horizontal black bars mark our two 
continuum regions. 
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Fig. 5. — Continuum polarization evolution in SN 2006ov. The 
black stars are the KAIT unfiltered light curve of SN 2006ov, with 
the solid line connecting the points to guide the eye. The blue 
diamonds and solid red circles represent the polarization in our 
two continuum windows. The end of the plateau, t v , corresponds 
to JD = 2454095. The plotted error bars on the polarizations 
include a 0.04% systematic error contribution. The two horizontal 
black bars mark our two continuum regions. 
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Fig. 6. — Spectropolarimetry of SN 2006my on day 38, plotted 
as 50 A bins. The Stokes parameters have been rotated by 174? 1, 
as described in the text. The solid line in the top panel is the total- 
flux spectrum. Note the high continuum polarization, near 1%, at 
continuum wavelengths near 7000 and 8000 A. The two horizontal 
black bars mark our two continuum regions. 

certainty of 0.04% in each Stokes parameter. Just as 
for SN 2006ov, the data for SN 2006my are consistent 
with negligible ISP in the host galaxy. In this object, 
the continuum polarization is distinctly lower than in 
SN 2006ov, although it is unclear how much of the dif- 
ference is due to the later epoch of the SN 2006my ob- 
servations. 

3.3. SN 2007aa 

Unlike the other two objects, SN 2007aa was discovered 
shortly after explosion, allowing us the opportunity to 
obtain data at early times and have a comparison for the 
late-time data. 

Our first epoch of spectropolarimetry was obtained 
on day —50, well before the end of the plateau, and 
only 28 days afte r disc overy. As discussed above, 
iFolatelli fc MorrelH (|2007| ) described a spectrum taken 
one day after discovery as being similar to that of 
SN 1999em at 20 days after explosion, implying an 
epoch for our first observation of ~ 47 days after ex- 
plosion. We used the S uperNova IDentification code of 
iBlondin fe TonTvl (f2007h to cross-correlate our day —50 



spectrum with a library of supernova templates, holding 
the redshift fixed. The top three individual SNe with 
spectral matches (ignoring spectra on multiple dates per 
object) were SNe 2004et, 2004dj, and 2006bp at 56, 56, 
and 58 days (respective l y) after their estimated explosion 
dates (jSanu et al.ll2006t)Vink6 et alJl2004 iQuimby et al.1 
l2007h . 

These estimates are all consistent with SN 2007aa hav- 
ing a normal plateau duration of ~ 100 days, with our 
first set of data being taken near the middle. Our second 
epoch, at day —22, was taken late in the plateau stage, 
but before the beginning of the steep decline phase, while 
our third epoch, on day 3, was taken during the steep 
decline and just before the beginning of the radioactive- 
decay tail. 

The two plateau epochs of spectropolarimetry are plot- 
ted versus wavelength in Figure [7J and in the q — u plane 
in Figure [H We have applied a rotation of 172° to the 
data in Figure [7] in anticipation of the late-time results. 
A low level of continuum polarization at early times in 
SN 2007aa is clear from Figure [7J At both epochs, mrsp 
was near zero and grtsp ranged over — 0.2% in the con- 
tinuum. Line polarization features are clearly present 
at the Balmer, Na I, Fe II, Ca II, and O I lines. The 
strongest line polarization feature, at Ha, had an extent 
of about 0.4% in the first epoch and grew to ~ 1% in the 
second epoch. Close inspection of ursp in the second- 
epoch data shows small variations at the positions of the 
line features. This is a hint that the line features indicate 
deviations from axisymmetry. 

To examine this in more detail, we have plotted the 
data points in the Ha line from day —22 in Figure [8] as 
the solid black circles. The data points are 50 A bins in 
wavelength from 6400 to 6750 A, representing a velocity 
range of —7500 to 8400 km s _1 relative to Ha. The 
radius of each circle is proportional to wavelength, so the 
smaller points are in the blue wing of the profile and the 
larger circles are from the red wing. The Ha points show 
a spread in both q and u in a manner inconsistent with a 
simple axisymmetric geometry. In axisymmetry, all the 
points would lie along a line in q — u space. Instead, the 
points show a q — u "loop." As one moves across the Ha 
profile from red to blue, one moves in a counterclockwise 
direction around the loop. 

Such loops were first ident i fied in supernova spectropo- 
larimetry by ICropper et al.l (|1988h in the Ha and Ca II 
near-infrared (NIR) triplet lines of SN 1987A. The in- 
terpretation of the loops in SN 1987A has remained ob- 
scure for over 20 years de spite the high quality of the data 
(jWang fc W heeler 2008) . Subsequent investigations have 
found q — u loops in systems as different from SN 1987A 
as SNe Iln (iHoffman et al.ll2008D . stripped- envelope SNe 
(iMaund et al.l l2007aBb|)7 and even SNe la fiC asen et al.l 
20031 IWang et alJl2003al IChornock fc Filippenkol 12008 



Kasen et al.1 ((2003) found that loops were a generic 



consequence of systems containing multiple components 
with misaligned axes of symmetry and constructed a few 
examples with very different geometries. The relevance 
of these models to SN II atmospheres during the plateau 
stage is unclear as both the photosphere and Ha-forming 
regi on are w i thin t he hydrogen envelope of the progeni- 
tor. IChugail ()1992l ) has a more promising model, specif- 
ically designed for SN 1987A, which interprets the Ha 
loops as evidence of a small number of individual clumps 
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Fig. 7. — Early-time spectropolarimetry of SN 2007aa obtained 
while the object was on the plateau, on days —50 and —22. The 
polarization data have been rotated by 172? 4 to align <7rsp with 
the direction of the late-time continuum polarization. The blue 
solid lines in the background of each 9fisp plot are the total-flux 
spectra from each date. Note the weak continuum polarization 
with significant line polarization features. The "droop" in mrsp to 
negative values in the 2007 March 18 polarization data at wave- 
lengths greater than 8200 A is unexplained and may not be real, 
although none of the standard or probe stars observed that night 
showed such an effect. In any case, none of the results presented 
in this paper rely on data from that portion of the spectrum. The 
horizontal black bar in each <Jr,sp panel marks our continuum re- 
gion. 

of 56 Ni creating pockets of high ionization and Ha exci- 
tation. At different radial velocities, we see the effects 
of different clumps causing the net polarization angle to 
rotate with wavelength. 

We defer additional analysis of the line polarization 
features in SN 2007aa to future work on the polarization 
of SNe II at early times because our focus here is on 
the continuum polarization and its temporal evolution. 
For now, we will note the following points. First, the 
long axis of the Ha loop on day —22 is mostly aligned 
with the P.A. of the late-time continuum polarization, 
possibly indicating a common symmetry axis. Second, 
the total polarization excursion in the line is ~ 1% at an 
epoch where the continuum polarization is likely to be 
very much lower, as discussed below. 

We will now discuss the continuum polarization in de- 
tail. We found our low-S/N late-time polarization data 
in the final epoch to be unreliable near 8100 A, so we have 
focused only on the bluer of our two continuum polariza- 
tion windows, 6800-7200 A. Figure [5] has our continuum 
polarization measurements plotted for each of our three 
epochs of observation. As was clear from Figure the 
observed continuum polarization on the plateau was very 
low. We measure formal polarizations of 0.11 ± 0.04% 
and 0.08 ± 0.05% on days —50 and —22, respectively. 



By contrast, the polarization after the end of the plateau 
jumped to 0.83 ± 0.12% on day 3. 

The late-time spectropolarimetry data of SN 2007aa 
from day 3 are plotted in Figure |H The data are ex- 
tremely noisy due to the faintness of the object at this 
time, but when binned up to 100 A per pixel the basic 
appearance of the data is qualitatively similar to that 
shown previously for SN 2006ov and SN 2006my. The 
continuum polarization is near 1% in <7rsp, but there 
are broad dips down to near zero polarization at the lo- 
cations of strong lines such as Ha, with ursp centered 
near zero. The polarization we measure at the peak of 
Ha is ( 9RSP , ursp) = (0.02%, 0.02%), with 0.1% error 
bars in each Stokes parameter, consistent with zero, as 
shown in Figure [8] 

A few of the dips in polarization at the lines appear to 
reach negative values of <?rsp , especially near 7650 A and 
8400 A, but these are only marginally statistically signif- 
icant. The wavelength bins that are the most negative 
are those with the least flux from the SN (due to absorp- 
tions from the telluric A band and the Ca II NIR triplet 
P-Cygni feature) and hence the lowest S/N, magnifying 
any systematic uncertainties in the sky subtraction. An 
integration over 8500 — 8800 A, representing the emission 
component of the Ca II NIR triplet, gives (<7rsp, ursp) 
= (-0.09 ±0.17%, -0.11 ±0.17%), which like the polar- 
ization at the peak of Ha is completely consistent with 
zero. 

If we interpret the late-time polarization of SN 2007aa 
as being similar to that of SN 2006ov and SN 2006my, 
then the contribution from ISP should be low in this 
object as well, given the low polarization in the lines. 
Further evidence for this can be seen in Figure 13 where 
we have plotted the peaks of the Ha lines for all three 
epochs (integrated over the same 6540-6600 A interval). 
The Ha peaks in the first and last epoch are near zero. 
The Ha polarization in the second epoch is higher, but 
still low (0.3%). That is also the epoch showing clear 
evidence of strong polarization throughout the line, so 
perhaps a non-zero polarization at the peak is to be ex- 
pected, partic ularly if an as pherical excitation distribu- 
tion is present (|Chugai|[T9 92). Small amounts of ISP have 
a similar effect to changing the origin in the q — u plane, 
so there is still a possibility of a pathological ISP cancel- 
ing out some intrinsic supernova polarization. However, 
taken together with the low continuum polarization mea- 
sured on the plateau, the simplest interpretation is that 
SN 2007aa, like SN 2006ov and SN 2006my, suffers min- 
imal ISP (< 0.1%). 

We have taken the continuum polarization measure- 
ments and plotted them versus time in Figure I10| along 
with the KAIT unfiltered light curve. The polariza- 
tion shows a strong increase from the plateau stage to 
the late-time observations taken during the rapid decline 
phase. Even if we are mistaken about the ISP being neg- 
ligible, the observation that there was a sharp change in 
the polarization associated with the end of the plateau 
is a robust result. This is strongly rem iniscent of the 
jump in polarization seen in SN 2004dj (|L eonard et al.l 
2006). In SN 2007aa, however, there was evidence of as- 
phericity in the form of line polarization on the plateau, 
while SN 2004dj did not even show any line polarization 
features before the sudden increase in continuum polar- 
ization. 
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Fig. 8. — Polarization data for SN 2007aa plotted in the q — u plane. The solid black circles connected by a line are from the day —22 data 
and are 50 A bins in wavelength from 6400 to 6750 A, representing a velocity range of —7500 to 8400 km s~ 1 relative to Ho. The points 
form a loop in the q — u plane. The radius of each circle is proportional to wavelength, so moving across the line profile from blue to red 
corresponds to moving around the loop in a counterclockwise direction. Also plotted are continuum polarization points (integrations over 
6800-7200 A) and a set of points representing the peak of the Ho emission line (integrations over 6540—6600 A). Note the low polarizations 
in the Ho peaks and in the continuum for the first two epochs, which together imply minimal ISP in this object. The day 3 continuum 
polarization represents a significant jump in polarization. The dashed lines represent the rotated 9rsp~ u RSP coordinate system, designed 
to have the gasp axis go through the day 3 continuum point. 



4. DISCUSSION 

The spectropolarimetry data presented above for the 
three objects after the end of the plateau are all qual- 
itatively similar with each other and with SN 2004dj 
((Leonard et al.ll2006h . The continuum polarization for 
each object is high, but flat and constant with wave- 
length between the strong line features in the red. This 
is a signature of the wavelength-independent nature of 
Thomson scattering. The lines produce large depolariza- 
tions and many overlapping lines significantly decrease 
the net polarization in the blue. The minimal polariza- 
tions at the peaks of Ha and lack of significant interstel- 
lar Na I absorptions are consistent with negligible ISP 
in all three objects, and therefore the large continuum 
polarizations are intrinsic to the SNe. 



Classically, the expected effect of lines forming in front 
of an aspherical electron-scattering photosphere is to pro- 
duce "inverse P-Cygni" line polarization features (|Jeffervl 
1989). Line emission near zero velocity produces a dip in 
polarization near the line center, and selective blocking 
of forward-scattered light produces a polarization peak 
in the blueshifted absorption trough. Higher polariza- 
tions can also be found in the red wing of the line as the 
opaque photosphere obstructs the observer from seeing 
unpolarized line emission from the far side of the ejecta. 
This simple picture provides an explanation for the line 
polarization features in photospheric-p hase spectrop o- 
larimetry of SN 1987 A and SN 1999em (|Jeffervlll991 allbl: 
iLeonard et alj |2001[) . and possibly our early data for 
SN 2007aa as well. However, we do not see emission 
peaks blueward of Ha or the Ca II NIR triplet in our 



10 



Chornock et al. 



11 A AT a. 

If" 


1 SN 2007aa 

A> A. /V f\K J \ ^ 




1 o ..." V 


fill J \ 

N PI ; 

rf ,l! 


Ra 







5000 6000 7000 8000 9000 

Rest Wavelength (A) 



Fig. 9. — Late-time spectropolarimetry of SN 2007aa from day 3, 
presented in 100 A bins. The solid line in the upper part of the top 
panel is the total-flux spectrum to guide the eye. High continuum 
polarization near 1% is seen in <?rsp with dips down near zero at the 
lines, while isrsp is consistent with zero. The spike in </fisp and dip 
in itRSP near 8100 A are spurious and are caused by systematically 
poor night-sky subtraction near that wavelength. The horizontal 
black bar marks our continuum region. 
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Fig. 10. — Time evolution of the continuum polarization of 
SN 2007aa compared to the KAIT unfiltered light curve (black 
stars connected by a line). The blue diamonds represent the con- 
tinuum polarization integrated over 6800-7200 A. The derived time 
of the end of the plateau, t p , is JD = 2454227.8. 

highest S/N late-time data (Fig. [3]). Instead, the polar- 
ization is simply low relative to the continuum through- 
out the line profile. There are a few contributing effects. 

The first is simple flux dilution. The continuum pho- 
tosphere and associated polarizing electrons are located 
relatively deep in the supernova ejecta and intrinsically 
unpolarized line emission from the optically thin outer 
layers of the ejecta travels unimpeded to the observer 
and dilutes the observed polarization from the contin- 
uum value. As an example, the [Ca II] A7300 doublet in 
the day SN 2006ov data represents emission from the 
low-density parts of the outer ejecta. The line emission 
strengthens over time as the ejecta thin out and the su- 
pernova begins the transition to the nebular phase. This 
line has a flux peak a little less than a factor of ^2 above 
the continuum (the continuum flux level is hard to de- 
fine at this epoch because of many overlapping P-Cygni 
features), and the polarization dips from 1.3 — 1.4% in 
the nearby continuum pixels down to ~ 0.8% in the line. 

The second reason is t hat photons scatter ing in lines 
are generally depolarized dHoflich et al. 1996). While the 
outer layers of the ejecta are optically thin in the contin- 
uum, it seems likely that the optical depth in the strong 



lines is still significant. The depolarizations over the full 
line profile from lines that are strong in the cool outer 
ejecta, such as Na I, Ca II, and Ha, are likely manifes- 
tations of this effect. Absorption and reemission in the 
line causes the photons to lose the directional informa- 
tion previously imprinted upon them by electron scat- 
tering. The net effect of a large number of overlapping 
lines in the blue is to r esult in a decrease in polarization 
at those wavelengths (jHowell et alJ 12001). In addition, 
the line-formation region is very geometrically extended 
compared to the electron-scattering photosphere, so the 
selective blocking of forward-scattered light is not an im- 
portant effect. 

A few of the strong depolarizations do not corre- 
spond to obvious strong features in the total-flux spec- 
trum. An example is the depolarization visible near 
5700 A in Figure El w hich was also present in SN 2004dj 
(|Leonard et al.ll2006f) . It is unclear which weak spectral 
feature is responsible for that depolarization and why the 
many other weak spectral features do not also result in 
similarly strong depolarizations. 

The time evolution of the polarization is of significant 
interest. SNe IIP generally do not exhibit strong polar - 
ization shortly after explosion ([Wang fc Wheeler! [2008). 
and yet SN 2004dj had a strong jum p in continuum polar - 
ization near the end of its plateau ([Leonard et al.ll2006h . 
The main motivation for this study was to determine 
whether SN 2004dj was unique, or if other SNe IIP ex- 
hibit similar polarization behavior. We have shown time- 
series plots for the two objects in our study with multi- 
epoch polarimetry in Figures [5] and [TQl but SN 2007aa 
was the only object with data taken on the plateau to 
establish the low early-time polarization. 

To make further progress, we searched for SN IIP po- 
larimetry in the literature. As discussed in §1, there 
are surprisingly few published datasets of SN IIP spec- 
tropolarimetry in the refereed literature. Since we are 
interested in the intrinsic continuum polarization, we 
must also exclude objects with large ISP corrections 
because estimates of ISP are un certain and require 
making assumptio ns. SN 1999em dLeonard et a.1.1 120011: 
Wang fc Wheeled l200l and SN 2004di dLeonard et all 
2006h are the only objects from the published literature 
to meet our criteria. Of the th ree objects in the study 
bv lLeonard fc Filippenkd ([20011) . SNe 1997ds and 1999gi 
both had large and uncertain ISP contamination, and 
SN 1998 A was a SN 1987A analog o f uncle ar relevance 
to our SNe IIP. lLeonard fc Filippenkol (|2005h also plotted 
continuum polarization versus time for several SNe IIP, 
but three of the five objects in their sample had ISP cor- 
rections of 1% or more. 

We collected photometry for the two objects from the 
literature and determined the date of the end of the 
plateau for each. The combined set of polarization versus 
time since t p is shown in Figure [IT] We must caution the 
reader that the data shown in this Figure are necessarily 
very heterogeneous. 

We also reanalyzed the published polarization data 
in order to make our comparisons as uniform as possi- 
ble. The published SN 1999em data generally do not 
extend sufficiently far to the red to include our favored 
continuum regions. We calculated the continuum polar- 
ization using as much of our 6800-7200 A window as 
was available in each observation from the datasct of 
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iLeonard et all (|200lD . The 1999 Nov. 5 and 2000 Apr. 
5/9 data included the full wavelength range, while for 
the 1999 Dec. 8 and 17 data we integrated from 6800 A 
to the red wavelength end of the data. The 1999 Dec. 17 
continuum polarization is the only one of the four val- 
ues to diff er from the synthesiz ed y-band polarizations 
quoted by ILeonard et al.l (|2001h by more than the error 
bars, possibly due to the limited wavelength range (the 
data from that night only extend redward to 6850 A). 
The data have not been corrected for ISP because the 
effect is believed to be small (~ 0.1%; Leonard et al. 
2001). 

As discussed above, applying the ve ry wide continuum 
region used by ILeonard et al.l (|2006l ) in their study of 
SN 2004dj to our SN 2006ov data would reduce the over- 
all polarization measured on day from 1.56% to 1.26%. 
To correct for this effect, we obtained the SN 2004dj data 
from D. Leonard and calculated the polarization in our 
two continuum windows in an identical manner to the 
other objects in this paper. Therefore, the peak con- 
tinuum polarization plotted in Figure [TT] is now 0.60%, 
~10% higher than th e published value of 0.56% given by 
ILeonard et all (|2006l) . 

Despite these caveats, we believe the basic trend 
shown in Figure [TT] is sound. At ear ly times, SNe IIP 
are s een to exhibit low polarization ((Wang fe Wheeleil 
2008). The large, extended, slowly rotating hydrogen 
envelopes of their progenitors are unlikely to be sub- 
stantially aspherical prior to explosion. Even an aspher- 
ical initial explosion will become more spherical when 
propagating outward through the remainder o f the star 
(|Chevalier fc Sokerl [l989t iKifonidis et all I2006T) . At the 
end of the plateau, the thick hydrogen envelope becomes 
mostly optically thin, allowing us to see deep into the 
ejecta and observe the aspherical core of the explosion. 
In addition, the observed polarization for a given geome- 
try is expected to be maximized when th e optical dept h 
to electron scattering is nearly unity (e.g., Hbflich 1991). 
To be fair, the dramatic jump in polarization was only 
directly seen in SN 2004dj and SN 2007aa. However, 
SN 2006ov does show a consistent increase in polariza- 
tion of 0.3% in both of our continuum regions over the 
course of just five days, indicating very rapid changes in 
the continuum polarization. 

It is instructive to compare to the supernova with 
the mo s t dense ly sampled polarimetric data, SN 1987A. 
Ueffervl (|1991bl ) presented ISP-corrected broad-band po- 
larimetry from a number of sources (jBarrettl 119881: 
iMendez et~ail 1 1 9881 ) ■ The BVRI polarizations do show 
a common trend of rising polarization at late times, with 
a possible spike in the B and V polarization percentages 
near the beginning of the radioactive-decay tail. While 
the SNe IIP in Figure [TTJ show declining polarization af- 
ter the start of the radioactive tail, the SN 1987A data 
continue rising up to ~ 1% at 20 0-300 days after ex- 
plosion. The Type lib SN 2001ig (jMaund et al.ll2007bft 
also showed a distinct jump in polarization associated 
with the hydrogen layer becoming optically thin and the 
photosphere receding into the He core. 

Prior to the publica t ion of the SN 2004dj polarization 
results, IChugai et al.1 (|2005| ) noticed unusual line pro- 
files present in that object in the nebular phase. The 
Ha, H/3, [O I], and [Ca II] lines all exhibited a double- 
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peaked structure with the main peak initially blueshifted 
by about 1600 km s _1 and moving redward with time, 
alon g with a simi l ar pea k on the red side of the pro- 
file. IChugai et al.l (|2005| ) modeled these line profiles in 
terms of a bipolar 56 Ni component protruding from a 
more spherical "cocoon " that produced an ionization and 
excitation asymmetry. IChugail (2006) found that mini- 
mal adjustments to the initial model could roughly match 
the polarization results at the second nebular epoch, al- 
though his models did not predict the decline in polar- 
ization after the beginning of the radioactive-decay tail 
seen in the data or the different polarization angle seen 
in the first nebular dataset relative to the later data. 

While asphcricities in the 56 Ni distribution are likely 
related to the polarization we observe in the three ob- 
jects in our sample, none of them exhibits late-time 
line profiles similar to those of SN 2004dj. The late- 
time Ha profiles of our objects are shown in Figure 1121 
SNe 2006my and 2007aa have mostly rounded Ha pro- 
files, with SN 2006my showing a small notch or dip near 
6575 A. The early-time spectra of SN 2006ov do show a 
pair of double notches, near 6535 A and 6590 A in the 
day —5 spectrum, corresponding to velocities of —1250 
and 1200 km s _1 with respect to Ha. The notches disap- 
peared in just five days and the line profiles after the end 
of the plateau were normal, by contrast with SN 2004dj, 
whose Ha line-profile asymmetries remained prominen t 
for at least 330 days after explosion (jChugai et alj [2005) . 
The blue and red notches in the Ha profile of SN 2006ov 
before the beginning of the radioactive-decay tail are 
reminiscent of the so- c alled "Bochum event" in SN 1987A 
(|Hanuschik fc Dachsl Il987f ). although at significantly 
lower velocities (~ 1200 km s _1 versus ~ 4000 km s _1 
in SN 1987A). Similar notches have been seen in other 
SNe IIP near the end of the plateau (e.g., at v w 
2500 km s" 1 in SN 1999em; Leonard et al. 2002a) and 
do likely i ndicate asphericities or clumps in the 56 Ni dis- 
tribu tion (|Chugail Il992t lUtrobin. Chugai. fc Andronoval 
fl995h . 

Given the evidence for a potentially clumpy 56 Ni distri- 
bution, it is worth considering whether clumps of 56 Ni or 
other deformations of the core-envelope boundary could 
be responsible for the continuum polarization that we 
meas ure, rather than an aspherical photosphere in the 
core. ILeonard et al.l (|2006f ) hypothesized that the change 
in the polarization angle seen from the first to the sec- 
ond epoch after the end of the plateau in SN 2004dj was 
due to some non-axisymmetric portion of the 56 Ni distri- 
bution beginning to be uncovered before the who l e core 
as the photosphere receded in the ejecta. IChugail (|2006l ) 
proposed an alternative explanation, which he termed a 
"spotty photosphere," invoking an inhomogeneous dis- 
tribution of the elements at the hydrogen-helium bound- 
ary due to fluid instabil i ties a t that interface during the 
explosion (Mii ller et al.l 11991] ) . These inhomogeneities 
would result in brightness and opacity variations that 
would be manifested as a nonzero net polarization that 
suddenly appears when the overlying layers become opti- 
cally thin and the photosphere passes through the core- 
envelope boundary, even if the core of the SN were es- 
sentially spherical on large scales. 

Our data for these three objects do not show any an- 
gle variations of the type seen in SN 2004dj. The polar- 
ization angles for SN 2006ov are identical for all three 
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Fig. 11. — Compi lation of SN I I P int rinsic continuum polarization measurements after correction for ISP. Th e SN 1999em data (purple 
squares) are from Leonard ct al. (2001) and the SN 2004dj data (cyan circles) are from Leonard et al. (2006), reanalyzed in a manner 
similar to the current objects (see text for details). The SNe 2006ov (blue stars), 2006my (orange diamond), and 2007aa (red triangles) 
data points are from this work. The absciss a represents time in d ays r elative to the end o f the plateau, t p , for each object, as defined in 
the text. We used i?-band ph otometry fromlLeonard et all (f200^l andlVinko et al.l ([200^ to derive t p for SN 2004dj of 2004 Oct. 12.4, 
and i?-band photometry from Leonard ct al. (2002a) and Hamuv ct al. (2001) were used to derive t p for SN 1999em of 2000 Feb 23.3. 
A systematic error contribution of 0.04% has been added in quadrature to all points to account for potential night-to-night instrumental 
variations, except for the SN 2004dj data, which already included a systematic error contribution. The black solid line is the well-sampled 
ij-band light curve of SN 1999em to guide the eye. The data are consistent with a sharp increase in the continuum polarization in the last 
10—20 days before the end of the plateau. 



epochs from day —5 through day 27. While the contin- 
uum polarization for SN 2007aa at early times is suffi- 
ciently small that the polarization angle is hard to re- 
liably measure (a small amount of ISP could result in 
dramatic changes in the inferred angle), the long axis of 
the q — u loop seen in the Ha line on day —22 (Fig. [8]) 
is approximately aligned with the late-time continuum 
polarization. The combination of low continuum polar- 
ization and significant line polarization (as seen in day 
—22 in SN 2007aa) may indicate that the Ha excitation 
is more sensitive to the underlying aspherical 56 Ni dis- 
tribution than the continuum and hence the alignment 
with the late-time polarization is not accidental. In addi- 
tion, the continuum polarization in SN 1999em remained 
at a constant P. A. over the course o f the first 161 days 
after explosion ((Leonard et alJ 12001). The polarization 
angle variation of SN 2004dj appears to be an outlier 
compared to these other objects, or else the first obser- 
vation after the end of the plateau was very fortunately 
timed to catch the polarization angle while it was rapidly 
varying. 

The simplest explanation for the constancy of the po- 
larization angles is that these SNe IIP have highly as- 
pherical cores surrounded by hydrogen envelopes that are 
more spherical. Thus, we conclude that our late-time 
continuum polarization measurements are most easily 



understood as good tracers of the underlying asphericity 
of the photosphere in the core, and not of a "spotty pho- 
tosphere" or other boundary effect that only produces a 
transient polarization signal when the photosphere enters 
the core. Observations taken well before the photosphere 
recedes through the boundary between the hydrogen-rich 
envelope and the helium-rich core already show signs of 
the asphericity that is present at late times. 

Our maximum continuum polarization of 1.56% in 
SN 2006ov on day implies a minimum axis ratio of 
1.45:1 in the co ntext of the as pherical electron-scattering 
atmospheres of iHoflichl (|1991h . if the object were viewed 
in the equatorial plane, and possibly higher if the sym- 
metry axis were less inclined to the line of sight. Despite 
the lower continuum polarization observed in SN 2004dj, 
iLeonard et al.l ()2006f ) derived a similar axis ratio for the 
core of that object after accoun ting for the preferred in- 
clination of IChugai etHI (|2005f ). 

SN 2006ov also shows some evidence for inner aspheric- 
ity in the late-time [O I] line profiles. The profiles of the 
AA6300, 6364 doublet in our latest spectrum (day 51) 
are shown in Figure HH Although the spectrum has not 
yet become fully nebular and the [O I] doublet does not 
yet completely dominate its portion of the spectrum, a 
double-peaked structure is apparent. The peaks in A6300 
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Fig. 12. — Late-time Ha line profiles. The spectra are labeled by 
their ages relative to t p . The bottom two spectra are of SN 2006my 
and 2007aa, while the remainde r are of S N 2006ov. The day —24 
spectrum of SN 2006ov is from ILi et all l l2007l) . while the day 51 
spectrum was acquired on 2007 Feb. 14.655 with LRIS using a 
setup similar to that for the other spectra. Two notches (marked 
with arrows) are seen in the earliest SN 2006ov spectra near 6535 
and 6590 A. At later times, none of the objects shows significant 
asymmetries in the Ha line profile. 

that this represents the oxygen distribution and not con- 
tamination by some other line can be seen from the pro- 
file of A6364, which shows a peak near zero velocity as 
well as one to the red; however, the redder peak in A6364 
is closer to 500 km s . 

Double-peaked [O I] profiles have been seen 
in stripped-envelope supernovae and attributed to 
a toroidal oxygen distribution (iMazzali et al.l 120051 : 
iMaeda et all 120081: iModiaz et all |2008|), but the veloc- 
ity splittings of the peaks in those objects are signif- 
ican tly higher. Alternative geometries are also possi- 
ble (jTaubenberger et al.|[2009t iMilisavlievic et al.ll2010D . 
Here, a more likely explanation is that we are seeing 
large-scale clumps or blobs in the inner oxygen dis- 
tribution, especially given that the main peak is at 
zero velocity. Such oxygen clu mps have been seen 
in st ri pped-envelope supernov a e (iFilippenko fc Sargentl 
119891 : iMatheson et all l2000MSilverman et al.l 12009) . 
SN 1987A (StathakisiiaDlilal, and the SN IIP 2002hh 
(|Pozzo et al.l I2006T K The clumps are usually attributed 
to Rayleigh- Taylor instabilities created at the interface of 
the oxygen-ric h and helium-rich l ayers in the core during 
the explosion (jMiiller et al.lll99"ll ). although other large- 
scale fluid motions during the explosion are possible. 

5. CONCLUSIONS 

We have presented late-time spectropolarimetry of 
three SNe IIP. All three objects had large continuum 
polarizations after the end of the photometric plateau 
(0-8 < P < 1.6%), while SN 2007aa had contrastingly 
smaller continuum polarization on the plateau (P w 
0.1%). These result s confirm that the po larization evolu- 
tion of SN 2004dj (jLeonard et al.ll2006l ) was not unique 
and that highly aspherical cor es exist in SNe IIP . This 
follows the trend identified by IWang et al.l (|2001h that 
the degree of polarization in core-collapse SNe generally 
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Fig. 13.— Late-time [O I] line profiles in SN 2006ov from day 51. 
The solid line shows velocities relative to A6300 and the dashed line 
shows velocities relative to the A6364 component of the doublet. 
Two peaks are observed in both, near and 850 km s^ 1 for A6300. 
These are likely a sign of large-scale clumps in the oxygen-rich 
inner core, perhaps as a result of Rayleigh- Taylor instabilities in 
the explosion. 

increases with time. 

It is important to note that an intrinsic continuum po- 
larization of 1.56% is quite high 5 . Alt hough there are ex - 
treme cases such as the SN Ic 1997X (jWang et al.ll2001l) . 
very few core-collapse SNe exhibit intrinsic continuum 
polarizations higher than this. Most stripped-envelope 
SNe have significantly lower continuum polarization (af- 
ter correction for ISP), including even extreme exam- 
ples such as the well-studied broad- lined S N Ic 2002ap 
ilKawabata et al.ll2002l: [Leonard et al.l2002UlWang et al.l 
2003b). The line polarization features in some objects 
are stronger, however. SNe IIP are the single most com- 
mon form of core-collapse supernova and represent the 
explosions of the most common massive stars, those at 
the low-mass end of the range of supernova progenitors. 
Our observations of these objects demonstrate that large 
degrees of asphericity are generic ingredients in even or- 
dinary core-collapse explosions and are not just limited 
to various unusual or extreme objects. 

This observational view parallels recent develop- 
ments in the theoretical modeling of core-collapse SNe. 
Modern simulations of supernova explosions have re- 
vealed a rich phenomenol o gy of multi-dimensio nal ef- 
fects (|Blondin et al.l 120031: Burrows et ail 120061) Ex- 
treme models invoking jets (jKhokhlov et al.lll999f) have 
had some success explaining the polarization evolution 
and other tracers of asphericity in co re-collapse SNe 
(jHoflich et al.l 120011: ICouch et all 120091) . although the 
mechanism by which jets would be launched in ordi- 
nary SNe IIP remains obscure, as jet formation is usu- 
ally thought to require atypically high angular momen- 
tum in the core (e.g., Burrows et al. 2007b). More- 
over, the current generation of neutrino-powered core- 
collapse simulations (e.g., Blondin & Mezzacappa 2007; 
Burrows et al. 2007a; Marek & Janka 2009) all generate 
complex structures in the explosion without needing to 
invoke jets. In particular, these simulations commonly 
create aspherical shocks dominated by low-order defor- 
mations that produce unipolar or bipolar outflows that 
seem promising for generating the core asphericities we 

5 We caution that the precise polarization value is dependent on 
the definition of the continuum regions. 
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infer from our polarization measurements. It remains to 
be seen whether these features persist in the upcoming 
more-detailed three-dimensional models. 

Unfortunately, most explosion models to date only sim- 
ulate the central core of the star and stop when an 
outward-going shock is launched (or has failed to launch). 
However, the shock continues to evolve as it traverses 
the star and our observations can only probe the final 
state of the core, which makes detailed comparisons of 
the observed asphericities to thos e seen in the simula- 
tions difficult. Only a few studies (|Kifonidis et al.l [20061 : 
iHammer et al.ll2009t ) have followed the aspherical shock 
wave as it propagates outward through the star. Future 
work in this direction, along with radiative-transfer cal- 
culations, is necessary to see if the existing models can 
reproduce the observations. 

It is interesting to note that SN 2006ov showed the 
largest continuum polarizations in our sample while 
it ejected the smallest amount of nickel and possi- 
bly h ad the lowest - mass progenitor. As mentioned 
above. ISmartt et al.l (|2009D quote a 56 Ni mass of 0.003 ± 
0.002 M©, a value about an order of magnitude lower 
than in most SNe IIP. That SN 2006ov ejected a low 
56 Ni mass can also be inferred by comparison of Fig- 
ures [5] and [101 SN 2006ov and SN 2007aa had similar 
magnitudes on the plat eau and SN 2006ov w as at a some- 
what smaller distance (|Smartt e t al. 2009), yet it was a 
magnitude fainter than SN 2007aa at the beginning of 
the radioactive decay tail. The progenitor mass limit for 
SN 2006ov is also ve ry low , M < 10 M© (jSmartt et al.l 
2009; ICrockett et al.l [2009). In addition, the expansion 
velocities of SN 2006ov are lower than those of the other 
SNe IIP (Figure [12]) . possibly indicating a low explosion 
energy. A connection between the low nickel mass, low 
progenitor mass, or low explosion energy and the high 
continuum polarization would be exciting if confirmed in 
other objects. 

We also encourage the enlargement of the published 
SN IIP spectropolarimetry sample, particularly of ob- 



jects lacking significant ISP. While observations at very 
late times are likely rare in the extant unpublished 
datasets, it would be fruitful to expand the number of 
objects plotted in Figure [TT] Only then will we be able 
to look for relationships between the core asphericity 
and other supernova properties. Trends with progeni- 
tor mass, explosion energy, and 56 Ni production will all 
be interesting to examine and may provide constraints 
on the theoretical uncertainties in the core-collapse su- 
pernova mechanism. 
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